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ABSTRACT 

We present near-infrared multi-object spectroscopy and JHKs imaging of the massive stellar content 
of the Galactic star-forming region W3 Main, obtained with LUCI at the Large Binocular Telescope. 
We confirm 15 OB stars in W3 Main and derive spectral types between 05V and B4V from their 
absorption line spectra. Three massive Young Stellar Objects are identified by their emission line 
spectra and near-infrared excess. The color-color diagram of the detected sources allows a detailed 
investigation of the slope of the near-infrared extinction law towards W3 Main. Analysis of the 
Hertzsprung Russell diagram suggests that the Nishiyama extinction law fits the stellar population of 
W3 Main best {E{J—H)/E{H—Ks) = 1.76 and i?K, = 1.44). From our spectrophotometric analysis of 
the massive stars and the nature of their surrounding Hii regions we derive the evolutionary sequence 
of W3 Main and we find evidence of an age spread of at least 2-3 Myr. While the most massive 
star (IRS2) is already evolved, indications for high-mass pre~main-sequence evolution is found for 
another star (IRS Nl), deeply embedded in an ultra compact Hii region, in line with the different 
evolutionary phases observed in the corresponding Hii regions. We derive a stellar mass of W3 Main 
of (4± 1) X 1O^M0, by extrapolating from the number of OB stars using a Kroupa IMF and correcting 
for our spectroscopic incompleteness. We have detected the photospheres of OB stars from the more 
evolved diffuse Hii region to the much younger UCHii regions, suggesting that these stars have finished 
their formation and cleared away their circumstcllar disks very fast. Only in the hyper-compact Hii 
region (IRS5), the early type stars seem to be still surrounded by circumstcllar material. 

Subject headings: HII regions - infrared: stars - stars: formation - stars: massive - techniques: 
spectroscopic 



1. INTRODUCTION 

Despite the impact on their surroundings, the forma- 
tion and early evolution of massive stars is poorly con- 
strained, primarily because of their scarcity and short 
lifetimes. OB stars are usually observed in very young 
star-forming regions, e.g. clusters and associations, 
which are still embedded. Near-infrared observations are, 
in most cases, the only way to study the stellar con- 
tent of these young regions. However, a pure photomet- 
ric characterization of young embedded stellar clusters 
is strongly hampered by highly varying extinction, un- 
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known distances and infrared excess of the young cluster 
members. 

On the other hand, spectroscopy offers an unambigu- 
ous identification of the massive stellar content and has 
proven to be a powerful method to find and character- 
ize the newborn OB stars. Stellar properties, like ef- 
fective temperature and luminosity, are derived based 
on the spectral features and extinction, distance and 
possible inf rared excess can be deter mined rel iably (e.g. 
Watson fc H anson 1997; Bl um et al.l[2 000: Ha nsoii et al.l 



2002t iB ik ct al .^2005i: iPuga et al.ll2006t iBik et al.ll20'Tor 

Puga e t al. 20 13). By comparing the effective tempera- 



ture and luminosity to stellar isochrones, the age as well 
as the mass of the stellar population in the young cluster 
can be derived m ore reliably than by photometry alone 
()Bik et al.l |2010[) . as the large extinction variations af- 
fects the identification of a clearly reddened cluster main 
sequence. Apart from the basic properties of the star 
clusters, also age spread inside the cluster, the evolu- 
tion of circumstcllar disks and the early evolution of the 
young OB stars can be characterized and studied. 

Analyses of increasing samples of stars in the same star 
cluster have produced evidence for age sprea ds in several 
cases. Based on the mas si ve sta r content, iClark et al.l 
(|2005f ): iNegueruela et al.l ()2010f ) find very little age 
spread (less than 1 Myr) in starburst cluster Wester- 
lund 1. However, in the Orion Nebula Cluster age spread 
has been found based on pre-main -scqucnce (PMS) stars 
()Da Rio et al.ll2010b[) . iKraus fc H illenbrand (2009) pre- 
sented evidence for an age spread in Taurus. Additional 



2 



Bik ct al. 



evidence for an age spread has been found in S255 where 
the most massive star is stiU deeply embedded and driv- 
ing an ou tflow, while the P M S population has an age of 
1-3 Myr (|Wang et aUlMll ). iKumar et all (pOOl found 
evidence for stellar clusters around high-mass protostcl- 
lar objects, suggesting a similar age spread. 

Additionally, young stellar clusters are the places to 
look for signatures of how massive stars have formed. 
Remnant accretion disks are being dispersed on a short 
timescale (jHollenbach et ahl 1 19941 ). However, massive 
Young Stellar Objects (YSOs) have been i dentified in 
several high-mass sta r forming regions (e.g. IBlum et al.l 
120041: iBik et al.l [20061 ). Studying the frequency of mas- 
sive YSOs as function of cluster age and mass will allow 
us to place time scales on the disk evaporation process. 

Before arriving on the main sequence, low- and inter- 
mediate mass stars arc still contracting and arc cooler 
and bigger than they will be on the main sequence. 
The intermediate-mass stars are detected as G and K- 
type stars and w ill evolve to A and F stars respectively 
()Bik et al.ll2010[) . Theoretical predictions show that stars 
more massive than 10 do not have an optically visible 
PMS phase (jPalla & Stabler 199(|). However, using deep 
near-infrared observations we can probe objects deeply 
embedded behind up to 50 magnitudes of visual extinc- 
tion and we might follow the PMS phase of the massive 
stars. 

Theoretical modeling of high-mass protostellar evolu- 
tion predicts that massive protostars might migrate to- 
wards the main sequence fo howing a path similar to in- 
termediate mass PMS stars (lYorke fc Bodenheimeiil2008l : 
iHosokawa fc Omukail [2009t iHosokawa et al.l 120 fofT Due 
to the high accretion rates, the massive protostars will 
swell up to a maximum radius of 125 Rq and exhibit 
a relatively cool temperature of ^5 000 K. In the HRD 
they will look li ke red supergiants (jMorino et al.l 119981 : 
iLinz et ahl 120091 ). After most of the stellar mass is ac- 
quired, the protostar is contracting to the main sequence 
and will gradually reach its main s equence size and tem- 
pcratm-e (Zinncck cr fc Yorkell2007f ). 

In this paper, we address these issues of massive star 
formation with near-infrared K-hand multi-object spec- 
troscopy as well as JHKg imaging of the massive stellar 
content of W3 Main, obtained with LUCI at the LBT. We 
perform for the first time a spectral classification of its 
massive stellar content, allowing the detailed assessment 
of the evolution ary status of these stars and their Hii re- 
gions. The W3 ()Westerhoutlll958| ) region is part of an ex- 
tended star formation complex in the Per seus spiral arm 
and is associated with W4 and W5 (see iMegeath et al.l 
120081 for an exhaustive review) , spanning an area of 200 
X 70 PC. W3 is located at a distance of 1.95 ± 0.04 kpc 
(|Xu et al.l[2006[ ). The W3 molecular cloud has a sheU 
morphology wrapped around two optically bright Hii re- 
gions; IC 1795 and NGC 896 of which IC 1795 has an 
age of 3-5 Myr ()Oevll2005l: iRoccataghata et al.ll20Tl[) . A 
much younger part of W3 is still deeply embedded in 
the molecular cloud surrounding the OB association and 
consists of 3 embedded clusters of compact Hii regions: 
W3-North, W3-0H and W3 Main. The latter region is 
the subject of this paper. 

Radio continuum observations reveal several com- 
pact a nd diffuse Hii regions in the W3 Main area 
(fWron-Willian^ [T97lL iHarris fc Wvnn- Williams! [19761) 



as well as ultra-compact Hii (UCHii) and hyper- 
compact Hii (HCHii ) regions (jClaussen et al.l 119941: 
iTicftrunk et al.l 
(fricftrun k et al.l 



199J. The l ocation of the C^^O 
T99I and NH3 (iTieftrunk et al.l [19981) 



molecular emission peaks in a radio quiet area south of 
the region W3 D. Sub-mm dust continuu m emission is 
found here too, as well as centered on IRS5 (jMoore et al.l 
[2OO7I) . 

Along with the first detection of the radio contin- 
uum sources, the first near-infrared imaging revealed sev- 
eral b right point sources associat ed with these Hii re- 
gions (I Wvnn- Williams et al.lll97^ . More sensitive imag- 
ing surveys (iHayward et al .ill 98a ITieftrunk et al.lll998l : 
lOiha et al.ll2004[ ) have revealed the stellar richness of the 
W3 Main region and identified all the candidate ionizing 
sources of the Hii regions via photometry. In the X-rays, 
where the extinction is much less severe, the entire PMS 
population is detected, and the cluster shows a spherical 
geometry with a size of 7 pc, about twice as large as seen 
in the near- infrared (jFeigelson fc TownslevI [20081 ). 

The different evolutionary phases of the Hii regions 
make W3 Main an ideal target to study age spread and 
the evolution of circumstellar disks around massive stars. 
The paper is organized as follows; in Section 2 we dis- 
cuss the LUCI imaging and spectroscopic observations 
and data reduction. Section 3 presents the near-infrared 
photometry, the spectroscopic classification and the de- 
rived HRD. In section 4 we discuss the evolutionary stage 
of the massive stars as well as the detected age spread, 
and summarize the conclusions in Section 5. 

2. OBSERVATIONS AND DATA REDUCTION 

Near-infrared observations of W3 Main were obtained 
with LUCI 1 (LBT NIR spectroscopic Utility with 
Camera and Integral-Field Unit for Extragalactic Re- 
search), mounted at the Gregorian focus of the Large 
Binoc ular Telescope (LBT) on Mount Graham , Ari- 
zona (iHih et al.l[2006l) . LUCI 1 ([Ageorges et al.l \20Wt. 
ISeifert et al.ll2010D~ is a near-infrared multimode instru- 
ment operating in seeing-limited imaging, long-slit spec- 
troscopy as w ell as multi-ob ject spectroscopy (MOS) 
mode (|Buschkamp et al.ll2010( ). 

2.1. Imaging observations 

W3 Main at a(2000) = 02''25'"37.5^ (5(2000) 
-t-62°05'24", was imaged in JHR^ with LUCI. The J 
and H data were taken on December 18 and December 
19, 2009, respectively, while the Kg data were obtained 
on Feb 5, 2010 with the N3.75 camera, providing a 4'x4' 
field of view at an image scale of 0.12"/pixel. 

The imaging data were taken with a detector integra- 
tion time (DIT) of 5 sec and a NDIT (number of integra- 
tions) of 10. A random dither pattern has been applied 
using 8 up to a maximum of 10 different positions and a 
jitter box of 30". The total exposure times per filter are 
listed in Table [1] Sky frames were obtained using the 
same instrument settings immediately afterwards, cen- 
tered on a off-cluster field at a (2000) = 02''25™41.0^ 
(5(2000) = +62°12'32". The final image quality of the of 
the J and H images is 0.6"and 0.9" for the Kg image. 

2.2. Spectroscopic observations 

W3 Main was observed with the LUCI MOS unit in the 
X-band on December 17 - 19, 2009. To ensure a proper 
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TABLE 1 

Observing details 



Observations 


Date (UT) 


DIT(s)^ 


Exp time (s) 


Objects 


J 


2009/12/18 


5 


800 




H 


2009/12/19 


5 


1000 






2009/02/05 


5 


900 




W3Kltl4K 


2009/12/19 


60 


1800 


IRS N3 


W3Kltl4Km2 


2009/12/18 


60 


1800 


IRS3a, IRS5, IRS Nl, IRS N6 


LSS 


2011/03/06 


60 


2400 


IRS4, IRS7 


WSBrightK 


2009/12/17 


30 


1080 


other objects 



^ Detector Integration Time 

alignment of the masks for the MOS observations, K- 
band pre- imaging was acquired on November 01, 2009. 
Three masks were created to cover the stellar content 
ranging from the most massive stars to the intermediate- 
mass PMS stars. We designed the masks such that the 
difference in magnitude between the stars in the mask 
was around 1 magnitude, ensuring comparable signal-to- 
noise ratio for all spectra in one mask. 

We used the 210_zJHK grating in combination with a 
slit width of 1.0", providing us with a resolution i?=4000. 

The spectra were collected with a spatial resolution of 
0.25"/pixel (N1.8 camera), providing the largest wave- 
length coverage (AA = 0.328 /xm). The MOS slits were 
typically 4" long, to allow for a 2" nodding for sky sub- 
traction while minimizing the contamination by other 
stars. 

In addition to the MOS spectra, we obtained a long 
slit spectrum of IRS4 and IRS7 on March 6, 2011 with a 
total integration time of 2400 sec. The instrument setup 
was the same as for the MOS spectroscopy. For the long- 
slit observations, the nodding offset was set to 45". The 
details of the observations are shown in Table [H which 
lists the different DITs and total exposure times used 
for the masks and long-slit spectra. The seeing of the 
spectroscopy data varied between 0.6"and 0.9". 

Telluric standard stars of early B type were observed 
immediately before and after the science frames. The 
standard star was observed using the same mask as the 
science observations, as well as being placed in two dif- 
ferent slits, covering the blue and the red part of the 
spectral band respectively. Immediately after the sci- 
ence and standard star observations, arc lamp and flat 
fields were taken to minimize the effects of flexure. 

2.3. Data reduction 
2.3.1. Imaging 

The near-infrared JHKg images were reduced with 
standard IRAP0 

routines, as described in lPasauali et ahl 
(|20Tl[) . The images where corrected for dark current 
and flat fielded using sky flats taken during morning 
and evening twilight. A sky frame was created by com- 
bining the images taken at the offset position using a 
weighted mean and rejecting the lowest 3 and highest 8 
values to remove the stars from the frames. The mas- 
ter sky frame was subtracted from the individual science 
frames. After background subtraction, the images were 

* IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 



corrected for geometric distortion. An astrometric so- 
lution was derived us ing stars in common with 2MASS 
(jSkrutskie et al.ir2006D . The corrected frames where com- 
bined and weighted by the average level of the back- 
ground, producing the final images (Fig. [1]). 

Photometry on the J HKg images was per formed using 
the starfinder software (jPiolaiti et al.|[2000l ). This pack- 
age is designed to fit an empirical point spread function 
(PSF) to a crowded steUar field. The stellar PSF was 
constructed using 22 bright and isolated sources in the 
frame. The radial extent of the PSF was set to 40 pixels, 
beyond the radius at which the PSF structure was dom- 
inated by background noise. The final PSF was created 
iteratively after selection and subtraction of neighboring 
sources from the image. Three iterations were performed 
before the final PSF was extracted. All sources at > 3 
a above the local background noise were included in the 
detected source list, and three PSF fitting iterations were 
performed to extract the faintest sources. 

We verified the results of our starfinder photometry by 
performing an additional run with the daophot package 
(jStetsoniri987f ) within the IRAF environment, daophot 
treats the background and the PSF model differently 
than starfinder^ allowing us thus to check our photomet- 
ric results from different approaches. Stellar sources were 
detected with the daofind task and aperture photometry 
was performed with the phot task in radii ~ 1 to 2 x the 
FWHM of the PSF. For each fiher a reference PSF was 
constructed by combining the PSF of at least 20 objects 
with the tasks pstselect and psf. PSF-fitting photome- 
try was performed with the allstar task, using the PSF 
models to fit all objects identied with a 3ct confidence 
level over the local background. Our results from both 
photometric runs are very similar, however, starfinder is 
able to better treat the saturated stars in the images. 

Several of the bright stars in the Hii region appeared 
saturated in our JHKg images. The J-band image has 
7 stars with counts above the linearity limit, while the 
H and Kg images have 8 stars in the non-linear regime 
of the detector. To obtain reliable magnitudes for these 
bright stars, starfinder was used to fit the wings of the 
saturated stars with the determined PSF to reconstruct 
the flux of the saturated part. However, a few of the 
brightest stars could not be recovered by starfinder due 
to their heavily saturated cores having negative counts. 
For those stars, the negative core of the PSF was removed 
and the peak of the PSF was flattened. The same area 
was also flattened in the reference PSF, and the best fit 
determined, such that the differences in the linear part 
of the star's and the reference PSF were minimized. 

This procedure is compared to the recovery done by 
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Fig. 1. — LUCI JHKs color image of W3 Main, North is up and East is left. The observed field of view is 4.5' X 4.5', corresponding 
to 2.55x2.55 pc at a distance of 1.95 kpc. The brightest region of extended emission in the north-east is the Hll region W3 A. The two 
diffuse, spherical Hll regions in the south are W3 K (south-east) and W3 J (south-west). In the center, the deeply embedded massive YSO 
IRS5 can be found. Annotated are the Hll regions Il l'ieftrunJt et al]|1997l ) as well as the candidate massive stars of which the spec t ra are 
presented in this paper (in italic). Additionally, the location of the 2 massive stars (NR386, NR559) identified bv INavarete et al.l l l20Tll ) 
are identified. 



starfinder on the non-saturated stars allowing the deter- 
mination of the uncertainties of the used method. The 
measured uncertainties where 0.1 mag for J, 0.07 mag for 
the i/-band and 0.08 mag for Ks- Comparison between 
the obtained magnitudes and the 2MASS values for the 
isolated stars shows a match within 0.1 mag for all three 
bands, suggesting that the above described method de- 
livers reliable values. 

Finally, we cross-matched the obtained catalogs for 
each filter to identify the sources that are detected in 
more than one band. We allowed a maximum off- 
set of 4 pixels (0.48") for the center of each point 
source. We calibrated the LUCI photometry with 
2MASS (jSkrutskie et al.ll2006D . The matching radius was 
4 pixels after transforming the 2MASS positions to the 
LUCI x,y coordinate system. Only stars with magni- 
tudes 10 < J < 16.5 mag, IQ < H < 14.5 mag, and 9.5 
< Ks < 14.0 mag were used to avoid saturation effects in 
LUCI and ensure reliable photometry in 2MASS at the 
faint end. As mismatches in data sets with different spa- 
tial resolutions and sensitivities are frequent, we selected 



stars with |mag(LUCI) - mag(2MASS)| < 0.5 mag from 
the median zero point of all stars in each filter. This led 
to a selection of 54 stars in J, 49 stars in H and 53 stars 
in Ks as final calibrators, yielding the zero points J^pt = 
29.205 ± 0.022 mag, i/^pt = 29.072 ± 0.014 mag, Ks zpt 
= 28.373 ± 0.018 mag. 

To obtain a realistic estimate of the photometric errors, 
we divided the datasets in two subsets with identical to- 
tal exposure time and created two auxiliary images for 
each filter. An empirical psf was constructed in the same 
way as on the deep dataset and starfinder photometry 
was derived for the auxiliary datasets. The photometric 
uncertainties are derived from the difference in magni- 
tude between the auxiliary datasets. For faint sources, 
not detected in the auxiliary datasets, the largest error in 
the same magnitude range, as derived for other sources, 
was applied as a conservative error estimate. 

The final photometric catalog, containing all the 
sources detected in the 3 bands, has a limiting magni- 
tude of '^21.5 mag in J, ~19.5 mag in H and ^^19.0 mag 
in Ks- A more detailed assessment of the completeness 
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limits will be done in a future paper where the entire 
stellar content of W3 Main will be discussed. 

2.3.2. MulU- object spectroscopy 

The MOS spectra were reduced with a modified version 
of lucired, which is a collection of IRAF routines devel- 
oped for the reduction of LUCI MOS spectra. First the 
raw frames were corrected for the tilt of the slit using a 
spectroscopic sieve mask as well as for distortion through 
an imaging pinhole mask. After correction, the slits were 
oriented horizontally on the frame and were extracted to 
be treated separately. The science data were corrected by 
a normalized flat field. The flat field was normalized in 
the spatial direction for every slit to remove illumination 
effects as well as the slit response. 

The long slit spectra were reduced using standard 
IRAF routines and no distortion or tilt correction was 
applied. The flat field was normalized by dividing with 
the mean of the flat fleld as the spatial variation in the 
counts of the flat fleld was very low. 

For both observing modes, the wavelength calibration 
was performed using the Ar and Ne arc wavelength cal- 
ibration frames. After the wavelength calibration, the 
sky backgrou nd was subtracted using the procedure by 
iDavied ()2007l ) to remove the OH line residuals. Finally, 
the ID spectra were extracted using the iraf task doslit 
and combined. 

The removal of the telluric absorption lines was per- 
formed in two steps. First the telluric absorption lines in 
the standard star spectrum were removed using a high- 
S/N telluric reference spectrurr0, providing a first rudi- 
mentary correction in order to reduce the absorption con- 
tamination of the telluric standard star. The Br7 absorp- 
tion line was fitted with a Lorentzian profile. The best 
fitting profile of Br7 was removed from the original stan- 
dard star spectrum, resulting in a atmospheric transmis- 
sion spectrum taken at the same airmass as the science 
observations. Finally, the IRAF task telluric was used to 
correct the science observations for the telluric absorp- 
tion lines. This procedure was applied to the standard 
star observations before and after the science observa- 
tions as well as in the two different slits. The correction 
showing the least telluric residuals was chosen as the final 
reduction. Finally, the spectra are normalized by fltting 
a spline function to the continuum. 

3. RESULTS 
3.1. Photometry 

The final reduced 3-color JHKs image is displayed 
in Fig. [H Apart from the stellar content, also several 
patches of extended emission can be seen. These areas 
coinci de with the location of the Hii regions in the radio 
maps (jTieftrunk et al.lll997f ). W3 A is the brightest Hii 
region in the near-infrared in the north-east of the im- 
age. The two diffuse Hii regions W3 K and W3 J can be 
identified in scattered light in the J-band due to the low 
extinction. Several dark dust lanes can be seen as well, 
showing that the extinction is highly variable. 

The Ks, H — Ks color-magnitude diagram (CMD) and 
J — H, H — Ks color-color diagram (CCD) confirm this 
large extinction variation over the observed fleld (Fig. 

^ obtained at NSO/Kitt Peak Observatory 



[2]). The CMD and CCD show all sources with a de- 
tection in all three bands and brighter than A's — 15 
mag. In the CMD, no clear reddened, cluster main se- 
quence can be detected. Apart from a foreground main 
sequence (at H — Kg w 0.2 mag), an enormous spread 
in H — Kg color of the sources in W3 Main is observed, 
caused by large extinction variations inside W3 Main. 
The CCD of the bright stars with Ks < 15 mag conflrms 
these huge extinction variations, as the location of the 
reddened main sequence stars (solid lines) is populated 
to at least Ak^ = 5 mag. Several sources are located 
to the right of the reddened main sequence, due to an 
infrared excess, most likely originating in a circumstellar 
disk (see Fig. [2]). Overplotted in the CMD is the loca- 
tion of the upp er end of the u n- reddened main sequence 1 
Myr isochrone (jLejeune fc Scl iaerer 2001) , with the spec- 
tral types annotat ed as well as the reddening vectors 
from the Jndcbetouw et al.l ()2005L hereafter Ind05) and 
iNishivama et al.l (|2009l . hereafter NI09) extinction laws. 
The candidate massive stars in W3 Main are expected 
to be located in between the reddening lines departing 
from a B3V and a 03V star respectively (diagonal lines 
i n the CMD) . 

lOjha et all (|2004[ ) identified 13 candidate OB stars 
which are bright near-infrared point sources without in- 
frared excess and associated with Hii regions detected in 
the radio (IRS2 - IRS N4 in Table [2]). The two candi- 
date massive stars (IRS4 and IRS7) have only H- and 
A's detections, but are added to the CMD, and are not 
plotted in the CCD. IRS5 is only detected in A'g. To 
identify massive Young Stellar Objects (YSOs) we added 
bright sources possessing an infrared excess and added 
two more bright sources in the fleld to complete the tar- 
get list for the spectroscopic observat ions (IRS N5 - IRS 
N8 in Table. [2]). Recently, Navaret e et al.l ()2011l ) pub- 
lished AT-band spectra of 4 sources in W3 Main. Two of 
them (NR386 and NR559) are added to our sample and 
included in our analysis. The two other stars, IRS N3 
and IRS N4. in common with our sample, are used as a 
cross check of the spectral type determination. 

A co mparison be t ween our LUCI photometry and 
those of lOjha et aD ()2004l ) and 2MASS shows signifl- 
cant differences. The LUCI photometry of the blue and 
iso lated sources is fu lly consistent with the photometry 
bv lOiha et all (|2004[) and 2MASS. However, for the red 
sources (IRS4 and IRS7) , as well as sources contaminated 
by strong background emission (IRS2a and IRS2b), a sig- 
niflcant difference (in H up to 0.7 mag and in Ks up to 
0.9 mag) is found. All these sources are loc ated inside 
Hii re gions and the lower resolution images of lOiha et al.l 
(|2004D as well as 2MASS possibly include more nebular 
contribution, altering the magnitudes. Additionally, for 
the extremely red sources, a possible small color term 
in the conversion between the different phot ometric sys- 
tems, for either our or the photometry of lOiha et ahl 
(|2004[ ). could also contribute to the observed differences. 
Due to the finer pixel scale and better seeing of our LUCI 
observations we consider our LUCI photometry more re- 
liable. 

The CMD (Fig. [2]) can be used to estimate the spec- 
tral type of the candidate massive stars by compari- 
son with the ab solute magnitudes predicte d by the 1 
Myr isochrone of iLeieune fc S chaerer* (20011) and adopt- 
ing a distance of 1.95 kpc (pCu et al..,2006. ). The 1 Myr 
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Fig. 2. — Left: H — Kb vs Ks color-magnitude diagram of W3 Main. Plotted as s quares are the candidate massive stars identified by 
lOiha et al.l II2004I ). in this work and the two stars identified bylNavarete et al.l II2011I ). Overplotted as a dashed line is the main sequence 
isochrone with an age of 1 Myr from lLeieune fc SchaereH I I2001I) with the location of the spectral types between 03 V and B3V are indicated. 
The two diagonal dotted lines show the redden ing lines for an 03 V and a B3V star using the extinction law of llndOSl . Additionally, the 
reddening vectors for Ak5=1 mag from llndOSI . (steep vector) and INI09I . (shallow vector) arc plotted. Right: J — H , H — Kb color-color 
diagram of W3 Main. The plotting symbols and lines are the same as in the CMD. The stars with an infrared excess, located more than 
3(T away from the reddened main sequence are plotted as triangles. 



TABLE 2 

Catalogue of the massive stars in W3 Main 



Object'' 



Hii a (J2000) <5 (J2000) 



region 



(h m s) 



J 

mag 



H 
mag 



Ks 
mag 



IRS2 


W3 A 


02 


25:44.3 


+62 


06:11.4 


12.23 


± 





10 


10.15 


± 





07 


8.97 




0.08 


IRS2a 


W3 A 


02 


25:43.3 


+62 


06:15.2 


12.90 


± 





03 


11.60 


± 





01 


10.65 


± 


0.02 


IRS2b 


W3 A 


02 


25:41.7 


+62 


06:24.2 


13.92 


± 





03 


11.71 


± 





01 


10.49 


± 


0.02 


IRS2c 


W3 A 


02 


25:47.1 


+62 


06:13.0 


13.58 


± 





03 


11.64 


± 





01 


10.59 


± 


0.02 


IRS 3a 


W3 B 


02 


25:37.8 


+62 


05:51.8 


19.55 


± 





03 


15.44 


± 





01 


12.56 


± 


0.02 


IRS4 


W3 C 


02 


25:31.0 


+62 


06:20.6 










16.45 


± 





02 


13.41 


± 


0.02 


IRS5 


W3 M 


02 


25:40.8 


+62 


05:52.3 


















13.16 


± 


0.02 


IRS7 


W3 F 


02 


25:40.5 


+62 


05:39.8 










17.07 


± 





02 


13.90 


± 


0.02 


IRS Nl 


W3 E 


02 


25:35.1 


+62 


05:34.5 


17.28 


± 





02 


13.91 


± 





01 


12.01 


± 


0.02 


IRS N2 


W3 H 


02 


25:32.6 


+62 


06:59.6 


11.74 


± 





10 


10.49 


± 





07 


9.88 


± 


0.08 


IRS N3 


W3 J 


02 


25:27.4 


+62 


03:43.2 


11.25 


± 





10 


10.35 


± 





07 


9.94 


± 


0.08 


IRS N4 


W3 K 


02 


25:44.8 


+62 


03:41.0 


10.31 


± 





10 


9.43 


± 





07 


9.01 


± 


0.08 


IRS NS'' 




02 


25:32.7 


+62 


05:08.1 


14.61 


± 





03 


12.75 


± 





01 


11.72 


± 


0.02 


IRS Ne'' 


W3 A? 


02 


25:47.4 


+62 


06:55.3 


15.04 


± 





02 


13.71 


± 





02 


12.97 


± 


0.02 


IRS 




02 


25:40.6 


+62 


05:46.8 


17.99 


± 





02 


13.30 


± 





01 


10.41 


± 


0.02 


IRS NS'' 




02 


25:26.3 


+62 


06:07.4 


20.99 


± 





30 


15.67 


± 





01 


11.91 


± 


0.02 


NR559<= 




02 


25:28.1 


+62 


05:39.6 


9.62 


± 





10 


8.53 


± 





07 


8.45 


± 


0.08 


NR386<= 




02 


25:37.5 


+62 


05:24.5 


12.36 


± 





03 


11.15 


± 





01 


10.39 


± 


0.02 



taken from lOiha et al.l l[200l ) 
^ naming con vention discussed in m ain text 
= Taken from lNavarete et al.l (20W) 

isochrone has been chosen as it likely represents the age 
of the region better than the ZAMS (see also Sect. 4 for 
a discussion on the age) . Using the relation b etween T^g 
and spec tral type from iMartins et al.l (|2005l ) for the O 
stars and lKenvon fc Hartmannl (jl995[) for the B stars we 
derived the matching spectral types. 

The so derived photometric spectral types are very 
sensitive to the chosen extinction law. Especially for 
the high-extinction sources, a different extinction law 
changes the spectral type significantly by up to 3-4 sub- 
types. We have tested several extinction laws and com- 
pared them to t he observed CCD ( Fig. [H . The ex- 
tinction laws of iRieke fc Lebofskvl ()1985[ ). llndOSi and 
iNIOQl provide a reasonable fit to the observed CCD, hav- 
ing a very similar E{J — H)/E{H — Kg) value of ~ 
1.73 (1.69, 1.73 and 1.76 , respectively). T he extinc- 
tion laws of lCardelh eFan (|1989h . lFitzpatrickl (|1999| ) and 



iRoman-Zuniga et al.l (|2007[ ) do not match the observed 
distribution of the stars in the CC D, showing a too sha l- 
low, E{J ^ H)/E{H - Ks) = 1.21 fCardeUi et al.l[T989l) 
and E( J- H)/E(H- K^) = 1.52 ()Roman-Zuiliga et all 
[20071) ■ or too stee p slope , E{ J - H)/E{H - K^) = 1.84 
(|FitzDatricklll999f ) of the reddened main sequence. In the 
CMD and C CD (Fig. ^ we have plotted the extinction 
law of llndOSl . 

The other parameter of the extinction law. the to- 
tal over selective extinction {R\ = A\/E{H — Kg)), 
is important when absolute values of Ak^ are derived 
e.g. to place the objects in the HRD or derive a pho- 
tometric spectral type. Rk is difficult to det e rmine 
from JHKg data only. IhdO^ . IRieke fc Lebofek^ (IT985h 
and INI09I provide different values for Rx, while hav- 
ing a similar slope of the near-infrared extinction law. 
Ihid05l derive Rk, = AkJE{H - K^) = 1.82 ± 0.09, 
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Fig. 3. — Normalized K-hand spectra of the massive stars in W3 Main as taken with the multi-object-mode of LUCI. Annotated with 
dashed lines are the spectral features which are important for classification of the stellar spectra. For the details of the classification see 
the main text. 



IRieke fc Lebofsk^ ([T985h find Rk, = 1.78, while iNIOl 
derive Rk^ = 1.44 ± 0.01. This parameter can only be 
constrained when the absolute magnitude of a reddened 
source is known. The effect of this parameter on the 
HRD will be discussed in Sect. 4.1. 

Table ini shows the result ing pho tom etric s pectral types 
for the extinction laws of llndOSi and INI091 For the low 
extinction sources, Ak^ < 1.5 mag, both laws result in 
similar photometric spectral types. For the highly em- 
bedded sources such as IRS4 and IRS7, the two spectral 
types differ dramatically, up to 4 subtypes. 

Apart from the candidate OB stars inside the Hii re- 
gions, other candidat e OB stars can be found in the 
CMD. iNavarete et all (|20llh identified one of the blue 



sources as a 07V star, not related to any Hii region. 
This would suggest that all the other bright stars could 
well be O stars, while the stars around /^s=13 mag and 
H — Kg ~ 2 mag are most likely candidate B stars and 
more widely distributed across the field of W3 Main. 

The CCD reveals a number of infrared excess sources. 
Taking into account the photometric errors, we found 15 
sources located more than 3a from the reddening line. 
The majority of the sources arc faint {Kg « 14 mag) and 
are most likely intermediate mass PMS stars showing 
circumstellar disk emission. Among the brighter sources, 
two are spectroscopic targets: IRSBa and IRS N8, with 
additionally IRS2a and IRS7 located very close to the 
3(7 line. The iiT-band spectra of these sources are thus 
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crucial to establish the true nature of their IR excess. 
3.2. Spectral classification 

K-hand spectroscopy of the candidate massive stars 
provides more reliable spectral types, based on the pres- 
ence of photospheric absorption and emission lines, than 
photometry alone. Additionally, it reveals the nature of 
the infrared excess sources detected in the CCD. 

The i^-band spectra of the massive stars are shown in 
Fig. m The majority of the spectra shows Br7 absorption 
with a narrow emission component. The absorption has 
a photospheric origin, while the narrow emission line is 
emitted by the surrounding Hii region. Due to variations 
of the nebular emission on small spatial scales, we made 
no attempt to correct the spectra for nebular emission. 
Other lines present in OB star atmospheres are He I 2.113 
/im, the N III complex at 2.115 /im (IRS2a and IRS3a) 
and He II at 2.1185 ^iu^ (IRS2). 

The LUCI spectra are compared with high-resolution 
reference s pectra of opti c ally c lassified O and e arly B 
stars from iHanson et all (|2005f ) and iBik et al.l (|20'05[ ). 
Reference spe ctra of mid- and late-B stars were added 
from iHanson et al.l (|1996[ ) . The reference spectra were 
degraded to the resolution of the LUCI spectra and ar- 
tificial noise was added to make the spectra comparable 
in S/N. A visual comparison is performed to determine 
the best matching spectral type. 

As an additional check, we determined the spectral 
types by comparing the measured equivalent widths 
(EW) of the Br7, He I and He I I absor ption lines with 
the observatio ns of IHanson et all (1 19961) and theoretical 
predictions of iLenorzer et al.l ()2004[ ). Since the Br7 ab- 
sorption lines are contaminated by the nebular emission, 
special care has been taken to measure the EW of Br7 re- 
liably. We fit a Moffat profile to the absorption wings of 
Br7 which are unaffected by the narrow nebular emission 
line, providing a more reliable EW determination. The 
Moffat fit resembles the observed line profile well and re- 
produces the line shapes of IRS N3 and IRS N5, which 
are unaffected by nebular emission . The two meth o ds de - 
liver very similar spectral types. iNavarete et al.l (j2011l ) 
obtained isT-band spectra of IRS N3 and IRS N4, and 
derived the same spectral types as we derive based on 
our LUCI spectra. 

As shown in Fig. [2l four stars with spectra (IRS2a, 
IRS3a, IRS N7 and IRS N8) show evidence for a near- 
infrared excess in the CCD. Two of them, IRS N7 and 
IRS N8, display an emission line spectrum. Additionally, 
the if -band spectrum of the deeply embedded source 
IRS5 also has these characteristics. These spectra show 
Br7 emission and IRS N7 also shows the CO bandhead 
emission at 2.3 /im, most likely originating f rom dense, 
hot circumstellar mate rial around the star (|Bik &: Thil 
l200l lBl um "eran[200l . These emission line spectra, as 
well as rising continua, suggest that these objects are 
massive YSOs, object surrounded by remnant accretion 
disk, possibly as a result of their formation ()Bik et al.l 
l200a . 

The other sources with a near-infrared excess, IRS2a 
and IRS3a, however, do not show a spectrum indica- 
tive of circumstellar matter. Their K-hand spectra show 
photospheric absorption and emission lines and they are 
classified as 08V-09V (IRS2a) and 05V-07V (IRS3a). 
Possibly, free-free emission of the bright Hii region could 



be the cause of their IR excess. This shows the impor- 
tance of spectroscopy to identify the stellar content reli- 
ably and disentangle the massive YSO candidates from 
the sources contaminated by nebular emission. 

Veiling by free- free emission makes the absorption lines 
more difficult to detect and could increase the uncer- 
tainty of the spectral type determination. In both IRS2a 
and IRS3a, the N III complex at 2.115 /im is detected, 
narro wing down the spect ral type determination to 05V- 
09V (jHanson et al.lll996l ) and making the effect of veil- 
ing on the spectral type determination less important. 

Using LUCI MOS spectroscopy, we have determined 
the spectral type of 13 OB stars in W3 Main as well as 
identified 3 ma. s sive Y SOs. Adding the two sources from 
INavarete et al.l ()201lD to our sample results in a total of 
15 OB stars. The final derived spectral types are given in 
Table[3]and vary from 05V - 07V for IRS3a to B3V-B4V 
for IRS N6. 

When we compare the spectral types with those de- 
rived from the photometry (Sect. 3.1), it turns out that 
the photometry usually predicts earlier spectral types. 
For some objects, like IRS2a, IRS3a, the spectral types 
derived with both methods agree well. However, in many 
cases the spectral types are rather different, e.g. for 
IRS2, IRS Nl and IRS N4. In the case of IRS N4, the 
photometric spectral type is 07V (for both extinction 
laws), for which the N III emission line at 2.115 /im and 
the He II absorption line at 2.185 /im are expected. How- 
ever, in the spectrum of IRS N4 these lines are not ob- 
served, resulting in a much later spectral type (B0.5V - 
B2V). 

3.3. HR-diagram 

After spectroscopic classification, the objects can 
be placed in the HRD. The intrinsic H — Kg col- 
ors a s well as the b o lomet ric corrections are taken 
from iMartins fc Pled (I2006D for the O stars and 
iKenvon fc HartmannI (11995 ) for the B stars. First, the 
E{H — Kg) is calculated and converted into Ak^. For the 
determination of Aks,, the two extinction laws ()Ind05l ; 
INI09D provide different values. In table [3J the values 
of Ak^ are given for th e Ilnd05l extinction law. To ob- 
tain the values using the lNI09[ extinction law, a factor of 
0.79x E{H - Ks) needs to be applied. 

After accounting for extinction, and applying the dis- 
tance modulus and the bolometric correction, the lumi- 
nosities of the massive stars in W3 Main are derived and 
they are placed in the HRD (Fig|4]). The left panel s hows 
the HRD with the extinction correction by Ilnd05l (i^K^ 
= 1.8 2), while the right diagram shows the correction by 
InIOQI (i?K, = 1-44). The error budget of the luminosity 
is dominated by the errors on the extinction coefficients, 
while the error in Toff is given by the uncertainty of the 
spectral classification, which in turn results in an extra 
uncertainty in the bolometric correction. 

Overplotted in the HRDs are the zero age main se- 
quence (ZAMS) and main sequence isochrones from 1 
to 3 Myr from iLeieune fc Schaered (|2001| ). the PMS 
isochrones for in termediate mass stars (M < 5Mq, 
iSiess et all l200Clf) as well as the theoretical birth- 
line (jPalla fc Stahleii Il990f ). Overplotted also are 
the high-mass protostellar evolutionary tracks from 
iHosokawa fc Omukail (|2009l ) for 10"^ and IQ-^Mq yr"i 
accretion rates. 
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TABLE 3 

Spectral types of the massive stars in W3 Main 



Object^' S/N ratio 


Photom. 


Spcctro. 






Class.'' 


Class. 


mag 



IRS2 


135 


04V/04.5V 


06.5V - 


07.5V 


2.32 


± 


0.39 


IRS2a 


125 


O9V/B0V 


08V 


- 09V 


1.90 


± 


0.34 


IRS2b 


130 


07V/09V 


BOV 


- BIV 


2.41 


± 


0.40 


IRS2C 


155 


O8V/B0V 


BOV 


- BIV 


2.08 


± 


0.36 


IRS 3a 


75 


05V/08V 


05V 


- 07V 


5.42 


± 


0.79 


IRS4 


100 


O6V/B0V 


08V- 


B0.5V 


5.69 


± 


0.83 


IRS5 


50 






YSO 








IRS7 


60 


O7V/B0V 


09V 


- B2V 


5.93 


± 


0.86 


IRS Nl 


105 


O8V/B0V 


B2V 


- B3V 


3.60 


± 


0.53 


IRS N2 


60 


09V/09V 


B0.5V 


- B2V 


1.28 


± 


0.27 


IRS N3 


110 


BOV/BOV 


BOV 


- B2V 


0.92 


± 


0.14 


IRS N4 


140 


07V/07V 


B0.5V 


- B2V 


0.94 


± 


0.24 


IRS N5 


80 


BIV/BIV 


BIV 


- B3V 


2.04 


± 


0.35 


IRS N6 


70 


B4V/B4V 


B3V 


- B4V 


1.47 


± 


0.22 


IRS N7 


110 


>03V 




YSO 








IRS N8 


70 


>03V/03V 




YSO 








NR559 




07V/07V 




07V 


0.33 


± 


0.20 


NR386 




O9V/B0V 


BOV 


- B2V 


1.54 


± 


0.23 



^ Spectral types derived with the llndOSi and IJN 109l extinction laws 

^ Va lues derived using the llndOSI extinction law, when adopting the 

IN109I law, the values need to be increased by a factor 0.79x E{H — K^) 




4.5 4.5 4.4 4.3 4.2 4.1 4.0 4.6 4.5 4.4 4.5 4.2 4.1 4.0 

log(T,fi) log(Teif) 



Fig . 4. — HRD of the massive stars in W3 Main. The dashed lines represent the main-sequenc e isochrones fromlLeieune Ik Schaereil 
l|20011 ) for 1, 2 and 3 Myr. The dotted lines show the PMS evolutionary tracks for a 4 and 5 Mq star l|Siess et al.ll2006l ) as well as high- mass 
protostellar evolution tracks for accretion rates of 10 ~* and 10 ~^ Mq yr~^ from [Hosokawa fc Omukai ( 200!^). The solid line represents 
the theoretical birth-line, above w hich a PMS star is not detectable in the optical i|Pana7r'5ta'h ler 1990). lef t: The stars are de-reddened 
using the extinction law of llndOSi extinction law. right: The same HRD, but using the extinction correction of lNIOgi . See also discussion in 
the text. 
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Comparison between the two HRDs shows that, on av- 
erage, the stars are clos er to the main sequence in the 
right diagram with the INI09I ex tinction correction. In 
the left diagram, with the llndOSi law, most of the sources 
are located to the right of the main sequence. Only 
IRS2a and IRS N6 fit better the main sequence in this 
diagram, while still in good agreement with the main 
sequence location in the lNI09l case. The very deeply em- 
bedded sources IRS4 and IRS7 are loca ted clo ser to the 
main sequence in the diagram with the INI09I law. Sev- 
eral sources, like IRS Nl, IRS N4 and IRS2 remain to 
the right of the main sequence in both diagrams. 

To clarify the description of the HRD, wc divide the 
HRD in three areas: the upper HRD with log L/Lq > 
5, where IRS2 and IRS3a arc located. The middle area 
of the HRD (3.5 < log L/Lq < 5), is where most of the 
stars are located. The lower HRD (log L/Lq < 3.5) is 
the area of the intermediate mass (PMS) stars where IRS 
N6 is located close to the PMS turn-on region. 

Two stars are located in the upper HRD, IRS2 and 
IRS3a, where the theoretical isochrones indicate that the 
most massive stars already show significant evolution af- 
ter a few Myr. As already noted in Sect. 3.1, IRS2 
has a significant difference between the photometric and 
spectroscopic spectral type. It is the most luminous ob- 
ject in W3 Main, but not the earliest type star (06.5V - 
07.5V). This is illustrated by its location to the right of 
the ZAMS, more consistent with the 2-3 Myr isochrones. 
The foreground extinction is not very extreme, therefore 
the location of IRS2 in the HRD docs not vary a lot for 
different extinction laws. 

IRS3a (05V - 07V), however, is very reddened 
(Ak,=5.42 ± 0.79 mag) and its location in the HRD 
depends strongly on the adopted extinction law. Addi- 
tionally, the presence of a near infrared excess (Fig. [2]) 
makes its luminosity more uncertain. However, the ef- 
fective temperature determination is more reliable and 
independent of extinction and excess. This allows us to 
place an upper limit on the age of 3 Myr, as stars on older 
isochrones will be of later spectral type than IRS3a. 

In the middle HRD, the position of the sources most 
consistent - within 3cr - with their expected main se- 
quence location. For spectral types between 09V and 
BIV, where most of the objects are found, the change in 
Teff and if -band magnitude is large between e.g. 09V 
and BOV. However, the the change in observable spectral 
features is not that dramatic. This is reflected in the 
large error bars of e.g. IRS N5 and IRS7. IRS Nl and 
IRS N4, however, are the furthest away from the main 
sequence. As their extinction is relatively low, their posi- 
tion docs not change much adopting different extinction 
laws. Below we discuss possible scenarios explaining this 
offset. 

A mis-classification of the spectral type is unlikely. 
Similarly as described in Sect. 3.2, the observed spectra 
of IRS Nl and IRS N4 are not compatible with the ex- 
pected spectral type derived from their luminosity (07V 
and BIV). 

Another explanation for the offset in the HRD is that 
the sources are at different distances. Calculating the 
spectrophotomctric distance of the stars, assuming that 
they all arc on the main sequence, would result in a 
distance between 1 and 3.4 kpc suggesting that subre- 
gions of W3 Main would be 1 kpc in the fore- or back- 



ground. iHever fc TerebevI ()1998D show that the entire 
W3/W4/W5 molecular cloud complex forms a connect- 
ing structure, both in spacial projection and velocity, 
suggesti ng that the entire com plex is at the same distance 
(see also lMegeath et aDl2008| ). Therefore, wc consider it 
highly unlikely that the offsets in the HRD can be ex- 
plained by different distances for the individual sources. 

Other effects, however, can be responsible for this offset 
in the HRD. Most of the massive stars are observed t o 
be binary stars (e.g. lBosch et al.ll200lHApai et al.ll2007| ). 
An equal mass binary would result in a mismatch in the 
observed luminosity by a factor of 2 (A log L/Lq = 0.3 
dex), with a multiple star system increasing the offset 
even more. Both sources would have to be a system of 
5-6 equal mass stars to explain their location above the 
main sequence, though. 

Stellar evolution can also be the reason why these ob- 
jects are located away from the main sequence. The loca- 
tion of IRS N4 would be compatible with a 10 Myr main 
sequence track, while IRS Nl would require an even older 
isochrone. This is highly unlikely, in particular for IRS 
Nl, as the star is located ins ide an UCHii region, which 
is still very young (- 10^ vrs. lWood fc"Chur chwcll 1989). 
IRS N4 is the central st ar of a diffuse and m ore evolved 
Hii region; as shown by iLada fc Ladal ()2003D a timespan 
of 10 Myr would be enough to disperse the surrounding 
gas. 

Another possible explanation for the HRD loci of the 
stars is that they are still in the PMS phase. Above 
the "bir th-line" no optically v isible PMS stars are ex- 
pected (|Palla fc Stahleii |1990[) . However, using deep 
near-infrared observations we might start to probe the 
PMS phase of massive stars. Theoretical modeling of the 
protostellar evolution predicts that high-mass protostars 
might have a similar evoluti on to intermediate PMS star s 
towards the main sequence (|Hosokawa fc Omukail l20091 . 

The exact path of the protostars in the HRD depends 
on the accretion rate. Overplotted in Fig. U are two 
tracks for accretion rates of 10~^Mq yr~^ and 10~^AIq 
yr~^. The location of IRS Nl would be consistent with 
an accretion rate between 10~^ and 1O~*M0 yr"-"^, while 
IRS N4 would require a slightly higher accretion rate. 

The location of IRS Nl and IRS N4 are consistent with 
the very last stages of the PMS evolution where the stars 
are contracting to their main-sequence temperature and 
luminosity. In this contraction phase, the accretion lumi- 
nosity only contributes 20% in the case of the 10"'' M© 
yr~^ and 10% 1O~^M0 yr~^ track respectively. The ab- 
sorption lines in the K-band spectra of IRS Nl and IRS 
N4 as well as the lack of near-infrared excess are con- 
sistent with the stellar contribution being the dominant 
source of luminosity. 

Surprising for those young objects however, is the lack 
of infrared emission coming from a circumstellar disk. 
The presence of a disk would make these objects similar 
to e.g. IRS N7 or IRS N8, which show an emission line 
spectrum and a near-infrared excess. An explanation of 
the lack of disk emission could be the fast dispersal of 
the disk by FUV photons of the central star or external 
FUV photons from other cluster members in W3 Main 
(see Section 4.1). 

Summarizing, the spectral classification of the mas- 
sive stars in W3 Main allows the comparison with stellar 
evolution models in the HRD. The location of the most 
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massive stars in the upper HRD suggests that IRS2 has 
an age of 2-3 Myr and IRSSa is slightly younger, while 
a speculative explanation of the displacement of IRS Nl 
and IRS N4 in the HRD could be high-mass PMS evo- 
lution. However, due to the large impact a change of 
extinction law has on the appearance of the HRD, it is 
difficult to assess with confidence that these stars are in 
the PMS phase. To prove this hypothesis, a good deter- 
mination of the stellar surface gravity, by means of mod- 
eling the absorption lines of high signal-to-noise spectra 
is needed to show that they have larger radii than main 
sequence stars. 

3.4. Energy budget of the Hll regions 

Apart from IRS N6, NR559 and NR386, all OB stars 
in W3 Main are located in Hii reg ions detected at radio 
wavelengths (jTieftrunk et al.lll997D . Assuming the radio 
emission to be optically thin, we recalculated the number 
of ionizing photons responsible for the radio flux using 
the distance of 1.95 k pc inst ead of 2.3 kpc quote d by 
Tieftrunk et"aI1 (fT997l) . With iMartins et al.l (|2005[) and 



Smith et all (I2002D . the number of ionizing photons are 



converted in a spectral type of a single star that would be 
responsible for the ionization of the Hii region (Table |4]). 
These spectral types are compared to the spectral types 
of the OB stars inside the Hii regions. In the case of 
W3 A (radio spectral type 06V), 5 OB stars are found 
in a single Hii region. Nevertheless, IRS2, the earliest 
type star, is the main source of ionization as the number 
of ionizing photons drops very fast with spectral type 
([Martins et al.ll2005| ). IRS2 is of spectral type 06.5V - 
07.5V emitting io48-44-48.80 Lyj-^aj^ continuum photons, 
while the 2nd most massive star, IRS2a (08V - 09V) 
only emits io4'i'-90~48-28 ionizing photons. 

For the Hii regions W3 A (06.5V), W3 B (07.5V), W3 
C (09.5V), W3 E (BIV) and W3 F (BOV), the most lu- 
minous stars provide sufficient Lyman continuum flux to 
ionize these Hii regions. For the diffuse Hii regions W3 H 
(09V), W3 J (08.5V) and W3 K (07.5V), on the other 
hand, the Lyman continuum flux of the brightest source 
is not sufficient to explain the radio flux. IRS N3 (BOV 
- B2V) in W3 J and IRS N4 (B0.5V - B2V) in W3 K 
are the brightest stars, located in the center of these Hii 
regions, and therefore are expected to be the main ion- 
izing sources. However, the expected spectral type from 
the radio flux is much earlier and inconsistent with the 
observed spectrum of the ionizing sources. Interestingly, 
the photometric spectral type of IRS N4 (07V) agrees 
better with that expected from the radio flux of W3 K. 
However, as discussed in Sect. 3.2, a 07V spectrum is 
not compatible with the the observed spectrum of IRS 
N4, so a miss classification of the spectrum is not likely. 
The spectroscopic cla ssification of IRS N3 and IRS N4 by 
iNavarete et al.l ()2011|) confirms our spectral classification 
based on a different data set. 

Leakage of FUV photons or dust attenuation would be 
able to explain w hy the ioniz ing star would be of earlier 
spectral type (e.g. iKurtz et al . 1994), but will not result 
in the fact that the observed radio emission cannot be 
accounted for by the ionizing star. iSmith et al.l (|2002D 
show that the number of ionizing photons is much higher 
(log Qo = 47.8 s~^) for a Bl super giant instead of log 
Qo — 46.5 s^^ for a Bl dwarf. This would come closer 



to the value derived from the radio emission. However, 
their position in the HRD (Fig |4|) does not comply with 
this explanation, since these stars are both too faint to 
be early B su pergiants, which wou ld have a luminosity 
of -lO'^-^i© (jMartins et"an 120051 ). while the observed 
luminosity is ~ (IRS N3) and - IO-^-^Lq for IRS 

N4. The presence of supergiants is also not consistent 
with the derived age of W3 Main. 



3.5. Cluster mass 

Determinations of the mass of embedded regions like 
W3 Main are hampered by the extreme extinction vari- 
ations observed in the cluster (Fig. [2]). No clearly red- 
dened main sequence is visible which cai i be compared 
to main sequence isochrone s directly (e.g. lAscenso et al.l 
120071 : iBrandner et al.ll2008() . However, the spectroscopy 
of the massive stars results in a direct mass estimate for 
the stars of the upper end of the mass function. 

The masses of the stars were calculated by compar- 
ing the stellar luminosities to the 1 Myr main sequence 
isochrone. Only for IRS2, the 3 Myr isochronc is used as 
it already evolved away from the location of the 1 Myr 
isochrone. This results in stellar masses of 35 Mq for 
IRS2 as the most massive star and 5 Mq for IRS N6, the 
lowest mass star in our sample. 

In order to obtain an extinction-limited sample, which 
is complete to the highest extincted massive star we have 
detected with spectroscopy, we took the mass of our red- 
dest OB star, IRS 7 (11 Mq) as our lower mass limit 
where we are complete. In this way, we select 13 out 
of 15 OB stars studied in this paper. Inspecting the 
CMD shows an additional 12 sources a bove this mass 
l imit without spectroscopic classification. INavarete et al.l 
(|2011|) identified one of the bright stars, not related to 
any Hii region, as a massive star, this suggests that the 
other stars might well be members of W3 Main. As- 
suming all these 12 sources are part of W3 Main, we 
end up with a total of 25 stars between 11 and 35 Mq. 
This number does not include the deeply embedded stars 
around IRS5 as they have even higher extinction. Addi- 
tionally, our observations only cover the central 2.5 pc x 
2.5 pc, while the cluster extends to 6 parsec b ased on X- 
ray observations (jFeigelson fc Townslevll2008l ). However, 
the Hii regions tracing the massive stars are all located 
in the central 2.5 pc. We calculate the mass of the clus- 
ter assuming that we are complete in the detection of the 
massive stars with < 5. 9 mag. 

Assuming a Kroupa IMF ()Kroupall2002f l with a slope 
r = -1.3 between 0.5 and 120 Mq and -0.3 for masses 
between 0.08 and 0.5 Mq and the total number of stars 
between 11 and 35 Mq (25 stars), we can extrapolate the 
mass function and obtain an estimate of the total stellar 
mass in the cluster. For the extrapolation we carried out 
Monte Carlo simulat ions assuming a random ly populated 
IMF as described in IBrandner et ahl (|2008f ). resulting in 
a total mass of (4 ± 1) x lO^Mg. 

The quoted uncertainty is only due to the effects of 
random sampling. When only the stars where spec- 
troscopy has confirmed their cluster membership are used 
(15 stars), then the derived mass is about a factor of two 
lower. Additionally, binaries are not taken into account 
as well as the deeply embedded population like IRS5. 
This would make the real mass likely higher. With 4000 
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TABLE 4 

Ionizing content of the Hii regions in W3 Main 



Hii 


log A^Ly 


Radio 


Stars 


NIR 


log Qo 


region 




Sp. type 




Sp. type 


s-l 


W3 A 


48.93 


06V 


IRS2 


06.5V - 07.5 


48.44 - 48.80 


W3 B 


48.43 


07.5V 


IRS3a 


05V - 07V 


48.63 - 49.26 


W3 C 


47.46 


09.5V 


IRS4 


08V - B0.5V 


47.00 - 48.29 


W3 E 


46.50 


BIV 


IRS Nl 


B2V - B3V 


< 46.10 


W3 F 


47.04 


BOV 


IRS7 


09V - B2V 


< 47.56 


W3 H 


47.85 


09V 


IRS N2 


BOV - B2V 


46.50 - 47.40 


W3 J 


48.16 


08.5V 


IRS N3 


BOV - B2V 


46.50 - 47.40 


W3 K 


48.45 


07.5V 


IRS N4 


B0.5V - B2V 


< 47.00 



Mq, W3 Main is about twic e as massive as the Orion 
Nebula Cluster (IO'Del]|l200lD . 

The total mass of the W3 GMC is estimated to be 
5x10^Mq based on ^^CO observations (jPickellQSdl ). but 
encompasses more than just W3 Main alone. Whe n look- 
ing at just the molecular gas towards W3 Main, iDickell 
(fl980f) derive a minimum mass for the high-density gas 
towards the W3 Main region of 5 x 10^M(7i, which is con- 
firmed by i^ciSQ observations (iThronsonI [l986D . To- 
gether with our mass estimate of 4 x IO^Mq, this would 
result in a estimate of the star-formation-efHciency (SFE) 
of ~44 %. Note that the southern part of the W3 Main 
region most likely already evacuated the gas, and there- 
fore this should be considered as a rough upper limit. 
Such high value of the SFE is consistent with the pre- 
diction of IBonnell et al.l (j2011[) for clusters forming in 
the gravitationally bound regions of a molecular cloud 
(SFE ^ 40 %). Wh ether the cluster will remain bound if 
the gas is removed (jPortegies Zwart et al.|[2010l ) remains 
uncertain and cannot be constrained with our current 
observations. 

4. DISCUSSION 
4.1. massive Young Stellar Objects 

Apart from the 15 OB stars, 3 massive YSOs have been 
detected. Their spectra are dominated by circumstellar 
emission and no spectral type can be derived from their 
i^-band spectra. However, based on their location and 
detection at other wavelengths, we can derive some upper 
limits to their luminosity. IRS5 is l ocated in a HCHii re- 
gion and is a small mult i ple system (IMegeath et al1l2005l : 
Ivan der Tak et al.ll2005t iRodon et al.ll2008[) with a total 
luminosi ty of 3 x lO^Lp) ioni zed by a handful of early 
B stars (jCampbell et alj|1995[ ). The two other massive 
YSOs in our sample, IRS N7 and IRS N8 are not as- 
sociated with a detectable Hii region. Also, they have 
no detectable sub-mm counter part. IRS N7 is located 
at the southern edge of the sub-mm clump centered on 
IRS5, wh ile IRS N8 is locate d east of the clump centered 
on IRS4 (IMoore et al. l [2007h . 

In contrast to the very embedded IRS5, both sources 
are detected in JHKs and therefore less embedded and 
probably more evolved than IRS5. The lack of radio 
emission would suggest that they are mid- or late B stars, 
as for more massive stars, an Hii region would be ex- 
pected. Similar o bjects have been det e cted in other star 
forming regions (iHanson et all [l997l: IBik fc Thil [200l 
iBlum et al.ri2004l: iBik et al.l l2006D ~but no clear reason 
has been identified for why these stars are still sur- 
rounded by circumstellar matter, while other OB stars 



already cleared their surroundings. IBik et ahl (|2006f ) 
showed that these objects have similar spectral and pho- 
tometric properties to Herbig AeBe stars, which are still 
surrounded by circumstellar disks even when they be- 
come optically visible. 

We detected OB stars in the diffuse, compact and 
UCHii regions. Only in the HCHii region W3 M, cir- 
cumstellar emission from IRS5 has been detected in the 
X-band spectrum. This shows that the accretion disks 
around the massive stars in W3 Main are destroyed 
within 10^ yrs. When an HCHii region evolves to the 
UCHii phase, circumstellar disks are destroyed by the 
strong FUV fi eld and stellar winds of t he recently formed 
massive stars (|Hollenbach et al-lfTOOl . 

4.2. Evolutionary sequence in W3 Main 

W3 Main harbors several different evolutionary stages 
of Hii regions, ranging from very young HCHii re- 
gions (few 10'^ yrs), U CHii regions (~ 10^ yrs 
IWood fc Churchwel]|[l98l to evolved, diffuse Hii regions 
(few 10^ yrs). All these regions are most likely formed 
out of the same molecular cloud. This provides the 
possibility to study the evolution of yo ung Hii regions 
and th eir stellar content in great detail. iTieftrunk et al.l 
(|1997f ) derived an evolutionary sequence for the Hii re- 
gions in W3 Main based on the morphology of the radio 
sources. The youngest are the HCHii regions W3 M and 
W3 Ca, with the UCHii regions W3 F, W3 C and W3 
E, slightly older, the compact Hii regions W3 B and W3 
A even more evolved, and the diffuse Hii regions, W3 K 
and W3 J being the oldest Hii regions in W3 Main. This 
classification can be compared to the ages of the massive 
stars deduced from their position in the HRD. 

We have detected OB stars in three diffuse Hii regions, 
two compact Hii regions and three UCHii regions with 
HCHii region W3 M harboring the high-mass protostar 
IRS5 (see Sect. 4.1). The position of IRS2 in the HRD 
suggests an age of 2-3 Myr, consistent with its location 
in a relatively evolved compact Hii region (W3 A) with 
a similar or younger age derived for IRS3a, harbored in 
a younger UCHii region. For the lower-mass stars (late 
O, early B) the isochrones are too close to each other 
to derive any age information. For IRS Nl, located in- 
side the UCHii region W3 E, the offset from the main 
sequ ence could be e xplained by high-mass PMS evolu- 
tion (jHosokawa fc Omukai 2009! ) . consistent with the ex- 
pected young age of the UCHii region. 

A similar sequence can be seen in the extinction to- 
wards the different Hii regions. The extinction varies 
from = 0.9 mag for the diffuse Hii regions W3 J and 
W3 K to very high extinction (Ak^ =5.9 mag) towards 
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W3 F. Only IRS3a, located in the compact Hii region 
W3 B does not follow the trend and shows a very high 
extinction (Ak^, = 5.4 mag). Such a sequence in extinc- 
tion indicates that the stars in diffuse Hii regions have 
ah'eady cleared out their surroundings and destroyed the 
molecular cloud, while those in UCHii regions are still 
embedded in their parental molecular cloud. 

Based on the presence of different evolutionary stages 
of Hii regions as well as the location of the most massive 
stars in the HRD, we can conclude that an age spread of 
a few Myr is most likely present for the massive stars in 
W3 Main. The presence of an older population of a few 
Myr is supported by the detection of several 10 s clas s 
III stars in the X-rays bv iFeigelson fc TownslevI ()2008D . 
These authors suggest that, at least some, of the OB 
stars have formed after the lower-mass PMS stars. A 
spectroscopic analysis of a sample of intermediate-mass 
PMS stars will provide more insights on this suggestion 
(Bik et al, in prep). 

A large number of young stellar clusters show evi- 
dence for an age spread, usually based on the analysis 
of the PMS population in the HRD. I n Orion an age 
spread of a few Myr has been found (jPalla fc Stabler! 
[T999HDaR"io et al.H2010bD . asweh as in the LMC in the 
young stellar cluster LH95 ()Da Rio et al.|[2010aD . In star- 
burst clusters, however, upper limits on the age spread 
of less than 1 My r have been found for Westerlund 1 
(jClark et al.l 120051: INegueruela et all 120101 ) based on the 
massive stars. Similar results are found for the PMS 
population of NGC 3603 and Westerlund 1 (Kudryavt- 
seva et al, in prep). This suggests that W3 Main, similar 
in mass to the Orion Nebula cluster, is not formed in 
one star formation burst as expected for starburst clus- 
ter, but, more likely, through a temporal sequence of star 
formation events. 

5. CONCLUSIONS 

We have obtained deep JHKs near-infrared imaging as 
well as ii'-band multi-object spectroscopy of the massive 
stellar content of W3 Main using LUCI at the LBT. We 
have confirmed the nature of 15 candidate OB stars by 
means of their absorption line spectra in the K-band. 
Three additional sources are identified as massive Young 
Stellar Objects based on their emission line spectra. 

From the analysis of the CCD and CMD, the extinc- 
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tion laws of Ilnd05l and INI09I provide the best fit of the 
slope of the reddened main sequence. Using either of the 
these two extinction laws the OB stars are placed in the 
HRD, with the iNIOOl providing a better match between 
the derived positions of the stars and the main sequence. 

Our most massive star, IRS2 (35 Mq) is already 
evolved away from the ZAMS and its location is con- 
sistent with an age of 2-3 Myr. Additionally, we find 
evidence for high-mass PMS evolution in the displace- 
ment of IRS Nl and IRS N4 from the main sequence. 
A total cluster mass of (4 ± 1) x 10'^ Mq has been de- 
rived by extrapolating the number of massive stars using 
a Kroupa IMF. 

The evolutionary sequence of the Hil regions seen in 
the radio morphology is consistent with a increase in 
extinction from the older (diffuse Hii) regions to the 
younger (UCHii) regions. We have detected the photo- 
spheres of OB stars from the more evolved diffuse Hii re- 
gion to the much younger UCHii regions, suggesting that 
the OB stars have finished their formation and cleared 
away their possible circumstellar disks very fast. Only 
in the HCHii phase (IRS 5), the massive star is still sur- 
rounded by circumstellar material. 

The above properties of the stellar content and Hii 
regions present evidence of an age spread of at least 2- 
3 Myr in W3 Main, similar to what has been observed 
in Orion and other young open clusters and in contrast 
to the single age population found in starburst clusters. 
This suggests a sequential star formation mode being at 
work for W3 Main, instead of a single star formation 
event. 
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